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ABSTRACT 

We discuss inflaton decays and reheating in no-scale Starobinsky-like models of inflation, calculat¬ 
ing the effective equation-of-state parameter, w, during the epoch of inflaton decay, the reheating 
temperature, T re h, and the number of inflationary e-folds, A*, comparing analytical approxima¬ 
tions with numerical calculations. We then illustrate these results with applications to models 
based on no-scale supergravity and motivated by generic string compactifications, including sce¬ 
narios where the inflaton is identified as an untwisted-sector matter field with direct Yukawa 
couplings to MSSM fields, and where the inflaton decays via gravitational-strength interactions. 
Finally, we use our results to discuss the constraints on these models imposed by present mea¬ 
surements of the scalar spectral index n s and the tensor-to-scalar perturbation ratio r, converting 
them into constraints on A*, the inflaton decay rate and other parameters of specific no-scale 
inflationary models. 
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1 Introduction 


A new generation of experiments on the cosmic microwave background (CMB), particularly 
Planck pj and experiments searching for 5-mode polarization, is providing detailed probes 
of models of cosmological inflation. In particular, recent data from Planck provide a very 
precise measurement of the scalar spectral index n s . Recent polarization results from the 
Planck and BICEP2 experiments [2] have focused attention on models [3] that predict 
relatively low values of the scalar-to-tensor perturbation ratio r, and a next generation 
of 5-mode polarization experiments is expected to produce results soon. Examples of 
low-r models include the Starobinsky model based on a R + R 2 extension of the Einstein 
Lagrangian OT, and related models such as Higgs inflation [?j, which typically predict 
v ~ 0.003. This is considerably below the current upper limit v 0.08 me). and models 
predicting values of r that are significantly larger than in the Starobinsky model may also 
be compatible with the data. 

We expect that the framework for physics at the Planck scale and below should be 
supersymmetric PGU. In addition to the myriad motivations from particle physics, su¬ 
persymmetry also renders technically natural the fact that the magnitude of the CMB per¬ 
turbations is small, by ensuring that radiative corrections to the requisite small mass scale 
and/or held coupling(s) are under control. The appropriate supersymmetric framework for 
cosmology is supergravity, but generic supergravity models of inflation soon encountered 
problems [12]. It was therefore proposed to consider models of inflation [1311151 based on 
no-scale supergravity m^\ B which are capable of mitigating these problems 0. 

Following the 2013 Planck data release, three of us re-examined [29b3Tj no-scale 
models of inflation based on a Kahler potential of the form 

I< = -3 In (t + T* - 

where T can be identified with the string compactihcation modulus and 0 is a generic 
matter held. With suitable choices of superpotential W(T, </>), no-scale models can repro¬ 
duce Starobinsky-like predictions with the inhaton identified as either the compactihcation 
modulus or a matter held, thanks to their conformal equivalence to R + R 2 gravity. There 
has subsequently been an outburst of interest in these and related no-scale models of infla¬ 
tion [32VI48] . In particular, we have shown how no-scale supergravity could accommodate 
models interpolating between the Starobinsky and chaotic quadratic models of inflation, an¬ 
alyzing their predictions for n s and r as functions of the number of e-folds during inflation, 
AT*, including also two-held effects [ 38l[39] . 

*We recall that compactifications of string theory lead generically to no-scale supergravity models fI5j . 
adding to their appeal. 

Wor some alternative supergravity-based models, see [TM25j . 
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We have recently studied various phenomenological aspects of such no-scale models of 
inflation, stressing how they could be embedded in compactifications of string theory US- 
We analyzed possible string assignments for the inflaton and matter fields, as well as mecha¬ 
nisms for supersymmetry breaking, inflaton couplings and decays. We showed that different 
no-scale supergravity models led to different estimates of the reheating temperature after 
inflation, T re h, and found a connection between the reheating temperature and the possible 
mechanism of supersymmetry breaking. 

The emerging data on n s and r are beginning to impose interesting constraints on 
the number of inflationary e-folds IV*, which depends on T re h and the equation of state 
during the epoch of inflaton decay, which is conveniently characterized by the effective 
equation-of-state parameter w int [451 .1491 150]. The cosmological data are therefore starting 
to impose supplementary constraints on inflationary models that may help discriminate 
among no-scale scenarios, casting some light on the mechanism of inflaton decay, and 
possibly supersymmetry breaking. 

In this paper we study these connections in some detail, comparing analytic and 
numerical calculations in Section 2 and evaluating w- mt , T re h and hence IV*- In Section 3 we 
apply these results in various Starobinsky-like no-scale inflationary models, and in Section 4 
we use the CMB bounds on n s and r to constrain IV* and thereby parameters in scenarios 
for no-scale inflation. Section 5 summarizes our results and discusses future prospects. We 
plan in a subsequent paper to study the low-energy constraints on supersymmetry breaking, 
and their complementary implications for no-scale models of inflation. 


2 On the Number of e-Folds in No-Scale Inflation 


In the slow-roll approximation and assuming entropy conservation after reheating, the 
number of e-folds to the end of inflation can be expressed as mmm 


IV* = 66.9 — In 


k * \ i 1 i 

^m) + 4 ln 


V? 


+ 


1 — 3 w 


int 


In 


MpP e nd / 12(1 T Wint) \Pend 


Preli 


12 hWreh i (2) 


where A;* is the wave number at the reference scale, ao and Hq are the present cosmological 
scale factor and Hubble expansion rate, respectively, 14 is the inflationary energy density 
at the reference scale, p en d and p re h are the energy densities at the end of inflation and 
after reheating, Wi nt is the e-fold average of the equation-of-state parameter during the 
thermalization epoch, and g ie h is the number of equivalent bosonic degrees of freedom after 
reheating: p reh = (7r 2 /30)sf re hT r 4 eh . 

We now discuss the evaluations of the quantities appearing in (J2|), with an initial focus 
on Starobinsky-like models of inflation that we extend later to related no-scale models. 
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2.1 The Inflationary Energy Density V* 


The Starobinsky potential V = |m 2 Mj,(l — exp V^" m p) 2 (where Mp = 1/\Z8 ttGn — 
2.4 x 10 18 GeV is the reduced Planck mass) is nearly scale-invariant for large values of 
the inflaton field 0: for 0 Mp, V — jm 2 Mp. This value is therefore a good first 
approximation to E- We can refine this value by recalling that the number of e-folds of 
inflation may be calculated in the slow-roll approximation as 


IV* 




1 r<f>e nd Y 

Mi 


76 

4M P 


(0end 


0 *) 


3 

4 



( 3 ) 

( 4 ) 


where 0* and 0 en d are the values of the inflaton field at the reference scale fc* and the end of 
inflation, respectively, and the prime denotes differentiation with respect to 0. Equation (141) 
may be inverted to obtain 0* in terms of the lower Lambert function W_\ (x). In practice, 
the asymptotic form W-\{x) = ln(— x) — ln(— ln(— x)) + • ■ ■ is sufficient to obtain a good 
estimate for 0*, namely 


~ \l -M P In 


-iV* 
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2 0end Jl+gg 
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This in turn implies that 


K ~ -m 2 M 2 P 


,— i /2~ ^end 

4iV* - + 3e v 3 m p 


( 5 ) 


( 6 ) 


In the range 50 < IV* < 70, this yields 0.728 m 2 Mp < E* < 0.734m 2 Mj,, a result that is in 
good agreement with the more exact values that we obtain from numerical integration of 
the equations of motion. 


The mass of the scalar field is not arbitrary, but is determined from the amplitude 
of the scalar power spectrum. At horizon crossing, the amplitude may be evaluated in the 
slow-roll approximation to be 

“ 127r 2 Mf0I/') 2 = 87T 2 (mp) Smh4 (vfefp) ' 


Using the approximation (]5]h this relation may be inverted to solve for the mass of the 
inflaton field, 


m ~ 8nMp 
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In the range 50 < IV* < 70 and using ln(10 ll M.5*) = 3.094 [1], this corresponds to 

1.218 < 10 5 (m/M P ) < 1.464. (9) 

Substitution in dHJ) leads to our final expression for the energy density at horizon crossing, 


K ~ 


I8n 2 As*Mp 
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( 10 ) 


which we use in our subsequent analysis. 


2.2 The Energy Density p e nd 

In the case of single-field inflation, the evolution of the homogeneous, canonically-normalized 
scalar 0 in the presence of a spatially-flat Friedmann-Robertson-Walker geometry is gov¬ 
erned by the equations 


0 + 3770 + U'(0) = 0 , 

^0 2 + E(0) = 3MpH' 2 , 


( 11 ) 

( 12 ) 


where H is to the Hubble parameter. Differentiating (1T2l) with respect to time and substi¬ 
tuting (flip yields the relation 


H = 


2 M 2 


(13) 


Using (gSD , the time dependence can be eliminated from the Friedmann equation, which 
leads to the Hamilton-Jacobi form of the equations of motion, 


[H’m 2 - FTTJW ) 2 = 


2Mp 


2Mp 


U(0), 


0 = -2M 2 P H\(p). 


The Hubble slow-roll parameters are defined by 


_o„2 W)V 


6ff(0) = 2 Mp ( tt/ 'A I — ei, 
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(14) 

(15) 

(16) 
(17) 


where ei t 2 are the first and second Hubble flow-functions, e± = —H/H 2 , e l+ \ = ii/(Hei) |{52| 
In terms of these parameters, the condition for inflation to occur is precisely 


a > 0 


< 1, 


(18) 
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which implies that inflation ends when e# = 1. 

Alternatively, one can consider the conventional potential slow-roll parameters 


en(0) = 


Ml 


W(0) = M 2 


V 

V'W 
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(19) 

( 20 ) 


which are fnlly determined by the shape of the inflationary potential. They can be expressed 
in terms of the slow-roll parameters via the relations 


W — £h 


3-77 H 


tH 


Vv = (2M|,ep) 1/2 —-f ( --— \ (e H + r ] H ), 

O — €h \o — e H ) 


( 21 ) 

( 22 ) 


which show that ey = 1 is only a first-order approximation at the end of inflation. It can 
be shown that the first term in ([22 ]) is of higher order in slow roll [51]. Neglecting this 
term, we can eliminate pu from equations (12T1 [22ll at the end of inflation, to obtain 


End of inflation: ey ~ (1 + \J 1 — r\y /2) 2 , 


(23) 


which can be used to calculate (j) en d- 

This equation involves the scalar potential and its first two derivatives, and can be 
solved in closed form in the case of a power-law potential V = a(0/Mp) n , yielding 

2n - 1 

2y/2 

This deviates from the exact result found by numerical integration of the equations of 
motion (1 11 Ml 2 p by less than 5% for n > 1. In the case of the Starobinsky potential, in a 
leading-order analytic approximation the end of inflation is reached when 

Starobinsky: 0 en d — In ^—(4 + 3\/3)^ Mp ~ 0.630Mp, (25) 

which is to be compared to the more exact value 0 e nd = 0.615Mp obtained by the numerical 
integration of the Friedmann and Klein-Gordon equations. 

The energy density at the end of inflation may then be obtained in a straightforward 
way by noting that the slow-roll parameter e# can be rewritten as ep = |(1 + w), where 
w = p/p is the equation-of-state parameter. When inflation ends, w = —1/3, which implies 

0e„d = V^end) ■ (26) 
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In the two cases discussed above, this may be evaluated to obtain 
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2 
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~ 0.182 m 2 Mp . 


( 27 ) 


The latter can be compared with p e nd = 0.175m 2 Mj,, which is obtained if we use the 
exact result for the Starobinsky potential, corresponding to the Hubble parameter H end = 
0.242 m. 


2.3 The Energy Density at Reheating p re h 

We calculate the energy density at reheating assuming that the inflaton decay is perturba¬ 
tive, with a rate T^,. As a first approximation, one can consider the decay to be complete 
when = t -1 . However, as we will see in Fig. |4j in general the decay of the inflaton is 
incomplete at this time. Instead, we assume here that reheating is complete when the bulk 
of the energy density is provided by the relativistic decay products of the inflaton: 

IT = Pl =1-5, (28) 

7 /V + /A 

for some suitable 5 <C 1. 

During reheating, the evolution of the inflaton field (ft and the relativistic decay prod¬ 
ucts can be described by the equations 


(ft + 3 H(ft + r 0 0 + V' — 0 , 

(29) 

P 7 477 p^y Y'fpPfp , 

(30) 

P<j> T P'y — 3MpH 2 . 

(31) 


It is only after integration of these equations that the moment when the decay is complete 
can be computed. However, we can find an approximate value when m T^ by averaging 
over the scalar field oscillations. The average energy density of the inflaton then corresponds 
to (p^) = (0 2 /2) + (V) — {(ft 2 ), and the average equation of motion (1291) simplifies to 

P<t> + 37/p 0 = —T^ . (32) 


These equations have the solution 


P(j) (f ) Pend 
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where u = T$tf, and we have assumed that the energy density of all relativistic degrees of 
freedom is negligible until the end of inflation [55]. If the decay rate is small, the reheating 
epoch spans a considerable interval of time, and t re g tend- In this limit, the scale factor 
and the Hubble parameter during the reheating epoch can be approximated as 


3 ( 1 +™) 


a(t) ~ a end ( — Pend (1 + w)t/Mp 


H ~ 


3(1 + w)t 


If we approximate w by its fzme-averaged value during reheating, 


(*^reh 


^eff — 


w(t ) dt , 


(35) 


(36) 


^end 


treh And 

then we can compute f re h by iteration. Consider as a first approximation w = 0, corre¬ 
sponding to the scalar held oscillations in the absence of decay. In this case, the solution 
can be combined with the constraint (I28|) to yield 


5 1 - 1 = e r ^ reh (r 0 t reh ) 2/3 




u 


2/3 e~ u du 


(37) 


= e r ^ reh (r^ reh ) 2/3 7 (|, r 0 t reh ), 

where in this context 7 denotes the lower incomplete gamma function. The relation (j37[) 
can be inverted numerically for any value of 6 < 1. For 6 < 10 -1 , the solution may be 
approximated by 

F^Iei- 0-754- 1.113 In5, (38) 

where the upper index denotes the degree of the approximation. An estimate for w e g may 
be derived by noting that the equation of state evaluated over the scalar held oscillations is 
just one third of the fraction of the total density provided by the relativistic decay products 
of the inhaton: 

i(0 2 )-(H) + i(p 7 ) p 7 /3 


(w) = 


= -H. 

P(j) P 7 3 


(39) 


\{<P) + (V) + (ft) 

Therefore, the zeroth order approximation to the time average of w can be calculated as 


(0) 

W eS ~ 


f r^reh 


7(i«) 


ST^tj-eh 


7 (|, u) + u 2 / 3 e~ u 


du ~ 0.271, 


(40) 


where for illustration purposes we have considered the end of reheating to occur when 
5 = 0.002 in ([28]) . The functional dependence of on 5 is illustrated in the left panel of 
Fig. [Ufor the case of an inhaton decay rate Y^ = 10~ 4 m. 

The computed value of w^g(5) may be substituted in (j35l) to calculate a hrst-order 
approximation to the reheating time, which may in turn be used to evaluate and 
so on. This iterative procedure relaxes after a few steps, resulting in 


rvS* =* 0.655 


1.082 In 5, 


(41) 











Figure 1: The effective equation-of-state parameters for = 10 -4 m as functions of the 
end of reheating defined by the parameter 5 in i\28\) . The numerical data are shown as red 
points, while the first-order and iterated approximations discussed in the text are displayed 
as dashed and solid lines, respectively. Left panel: The time average w e g. Right panel: The 
e-fold average wi nt . 


and (0.002) ~ 0.273. We have checked that the iterative solution is not sensitive to 
the initial choice of w. We see in the left panel of Fig. [Tjthat , shown as the solid black 
line, agrees very well with the results of integrating numerically the evolution equations, 
shown as the red points. Numerical results for w e s are shown in Fig. [2] as functions of the 
decay rate of the inflaton. We see that the results converge to the value w e s ~ 0.271 for 
<C m, a result consistent with the approximation (l40]h 

The energy density at the end of reheating may then be approximated by 


Pre h = 3 M 2 H 2 ch ~ 3 Ml ( , 2 u ) =U 1 + w eS )~ 2 M 2 P T 2 J0.655 - 1.082 lnh)" 2 . 

\3(1 + WefffirehJ 3 V 

(42) 

The corresponding reheating temperature T re h (assuming rapid thermalization [5£]) is given 
by 


F re h 


/ 30p reh \ 1/4 

V 7r2 .9reh ) 


(43) 


where the number of degrees of freedom g re h would be 915/4 for T re h above all the sparticle 
masses rh and falling at lower T re h, e.g., to g ie h = 427/4 for m t < T re h < fh. 
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Figure 2: The effective equation- of-state parameter w e s as a function of the decay rate of 
the inflaton. The data converge to the value w e s = 0.271 for m. 

2.4 Result for iV* 

Using the previous results, we can rewrite (J2J) for Starobinsky-like models in the form 


IV* = 68.66 - In (—+ \ In (4 5 *) - ^ In ( TV* - \/ -krF + ve v 3 m p 


\aoHo J 4 

1 - 3Wi n t 
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12(1 + Wi n t) 


(2.030 + 2 In (T^/m) - 2 ln(l + w eS ) - 2 ln(0.655 - 1.082 In 5)) , (45) 


\r /3 0end . 3 n c« 

A/ 2 = — \/ —|- _ e v 3 m p ~ 0.86 , 

M p 4 


(46) 


then (gap can be inverted in terms of the (upper) Lambert function Wq, resulting in 


JV. = -W„ (2e 2(M+X2 >) - jV 2 

— Ml + — In 2 — — In (2(Mi + Ml) + In 2) + • ■ ■ 


(47) 

(48) 


which is the basis for our subsequent analysis. 
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This expression for iV* depends explicitly on the decay rate of the inflaton, T^, and 
also implicitly, since its value affects the effective equation of state during thermalization, as 
characterized by the e-fold average parameter lUintQ-V) as we ^ as the ^ me average tc e ff(T 0 ) 
introduced previously. In the previous section we have derived an estimate for u> eff , finding 
that it has a universal value for m. The e-fold average of the equation of state u> int 

is given by the time average of w weighted by the Hubble parameter, 


-Ar reh (5) 


Wmt = 


w{ri) dn = 


^reh (^) 


Areh(^) ^end JN end N Ie h(5) ^end Jt end 

Following the same procedure for tc e fr, we can approximate U 7 i nt as 


( 0 ) 

7 /;> ' ~ 
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7(5, w) 


du 


w(t)H(t ) dt. 


0.731 


3 In' /3 


4Pend treh/Mp J 0 

0.743 


7 (|, u) + u 2 / 3 e~ u u ln(2.67m/T^) ’ 


mt ln(3.40m/r^) ’ 
for T^ m and 5 = 0.002. 


(49) 


(50) 

(51) 


The previous semi-analytical results can be compared with the results of numerical 
integration of the equations (129]) - (13TT) . The dependence on the parameter S of the effective e- 
fold -averaged equation-of-state parameter w int is shown in the right panel of Fig. Q] Iteration 
from the first-order analytic approximation to the e-fold-averaged parameter tCi nt does not 
converge as rapidly as that for the time-averaged parameter w e s (shown in the left panel 
of Fig. [T]and in Fig. [2]). We see that the 5-dependence of the iterated approximation w°f t 
for Y^/m = 10~ 4 (solid line) mirrors that of the numerical solution (red dots), though with 
a fractional offset < 10 %. 


Fig. [3] shows numerical values of w- m t together with the estimate (j5B as a function of 
Y^/m. Also displayed is a fit to the data, given by the equation 

0-782 

Wint ln(2.O96m/r 0 ) ' (5 } 

However, as already seen in the right panel of Fig. [T] the fractional difference between the 
numerical result and the iterated analytic approximation decreases as Y^/m —> 0 , as seen 
in the insert in Fig. EJ and is < 10% for the range of Y^/m of interest for our subsequent 
analysis. 


Fig. [4] illustrates the relationship between w(t) and ( w ), which corresponds, by virtue 
of (l39]h to 1/3 of the energy density in radiation. As noted earlier, we see explicitly that 
the estimate t re h ~ T / 1 (shown by the vertical green line) does not account fully for the 
decay of the inflaton into the relativistic degrees of freedom^. 


■'Note that we have chosen a large value of T^/m to be able to see graphically the oscillations as a 
function of mt. For smaller T^/m, the frequency of oscillations would be larger and the details of the 

oscillations would disappear. 
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Figure 3: The effective equation-of-state parameter w int as a function of the decay rate of 
the inflaton. The solid line corresponds to the fit K5E) to the data on the average of w, and 
the dashed line represents the estimate The inset displays the fractional difference 

between this approximate expression and the fit f5E) , which is small for -C m. 

0.6 

0.4 

£ 0.2 

0.0 

- 0.2 

mt 

Figure 4: The equation-of-state parameter w(t) for a inflationary model with Starobinsky 
potential and decay rate T^/m = 2 x 10~ 2 . The inflaton field oscillations are shown as a 
blue line, and the smooth (orange) curve interpolating between zero and 1/3 corresponds to 
the running average density ratio ( w ) f3fh) . The vertical (green) line is located at the point 



Fig. [5] shows the time at which reheating ends as a function of the decay rate. We 
see reasonable agreement between our analytical approximation (1411) for 5 = 0.002 (solid 
line) and exact results found by numerical evaluation of the equations of motion (l29]h(l3Tjh 
which are represented by blue dots. The energy density at the end of reheating is displayed 
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in Fig. [6l together with the approximation (1421) with the value of w e & that is determined 
by numerical integration. It is evident that the approximate expression is a good fit for the 
data, with a deviation < 3% for ras shown in the insert in Fig. [6l 



5 x 10 -5 1 x 10 -4 5 x 10 -4 0.001 

r 0 /m 


Figure 5: A comparison between numerical and approximate analytical calculations of the 
end of reheating as a function of the inflaton decay rate The (yellow) dots are obtained 
from the numerical integration of the equations of motion f2fA) - f3l\) . The estimate ffl\ ) 
with 5 = 0.002 is represented by the solid line. The inset displays the fractional difference 
between the approximate expression t[ff( and the exact numerical result. 


3 The Number of e-Folds in Representative No-Scale 
Inflation Models 

The preceding Section shows that we have good numerical and analytic control over the 
inflaton decay and reheating process, which we now use to calculate the number of e-folds 
IV* in some representative no-scale models of inflation. 

We see from (1441) that IV* depends on both explicitly and implicitly via the de¬ 
pendences in w int and w e s, which have been shown in (|52|) and Figs. [2]and|3j We use these 
in the general expression (144|) to calculate IV* as a function of For this purpose, we 
use the Planck pivot point fc* = 0.05/Mpc, corresponding to k^/aoHo = 221, and take the 
MSSM value of g re h = 915/4. Fig. [T] displays the calculated value of IV* over a wide range 
of parametrized by 

lyl 2 

r 0 = mf-, (53) 
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rym 

Figure 6 : Comparison between numerical and approximate analytical calculations of the 
energy density at the end of reheating as a function of the inflaton decay rate The 
(yellow) dots are obtained by numerical integration of the equations of motion (HPj)-(EI|). 
The solid line corresponds to the approximation with the value of w e s found by nu¬ 
merical integration. The inset displays the difference between the approximate expression 
with w e s 7 ^ 0 and the exact numerical result. 

with a coupling ranging from y — 1 to the value y ~ 10 -16 , in which the latter would corre¬ 
spond to a reheating temperature T re h — 10 MeV, below which the successful conventional 
Big Bang nucleosynthesis calculations would need to be modified substantially. Within this 
overall range, we discuss the values of TV* found in specific no-scale models whose inflaton 
decays were discussed in [30] , 

3.1 Decays via Superpotential Couplings 

In one class of model discussed in [40], the inflaton was identified as an untwisted matter 
field in some suitable string compactihcation, with direct decays to matter particles via 
a Yukawa-like superpotential coupling. One possible realization of this scenario would be 
further to identify the inflaton as a singlet (right-handed) sneutrino N with a couplings 
y u HLN to light Higgs and lepton doublets [31j[57]. The perturbative decay rate of such 
a sneutrino inflaton would be given by (1531) with y identified as the neutrino Yukawa 
coupling y v . We use this as a representative of the broader class of matter inflatons that 
decay directly to matter particles via trilinear superpotential couplings. 

Within the sneutrino inflation scenario, one might wish to consider values of y u < 1, 
the upper limit corresponding to a value of the Yukawa coupling similar to that of the top 
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Figure 7: The values of N* in no-scale Starobinsky-like models as a function ofT^/m, 
for a wide range of decay rates. The diagonal red line segment shows the full numerical 
results at 5 = 0.002 over a restricted range of T^/m, which are shown in more detail 
in the insert, and the diagonal blue strip represents the analytical approximation (O for 
io- ' 3 < 5 < 10 1 . The difference between the results from evaluating w- mt via the iterative 
procedure and through the analytical approximation in the fit T52 1) are indistinguishable in 
the main plot, but are visible in the insert, where the solid black line corresponds to K 5H) . 
The right vertical axis shows the values of n s in Starobinsky-like no-scale models, for which 
the tensor-to-scalar ratio varies over the range 0.0034 < r < 0.0057 for IV* in the displayed 
range. The vertical coloured lines correspond to the specific models discussed in Section 3, 
and the horizontal yellow (blue) lines show the 68 and 95% CL lower limits from the Planck 
2015 data, which vary slightly in related no-scale models, as discussed in Section f. 
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quark. In this case, we estimate 


y — 1 : IV* ~ 55.5, (54) 

as shown by the vertical red line in Fig. [71 However, such a large value of y would reheat 
the Universe to a very high temperature T reh ~ 10 14 GeV, which would lead to an overpro¬ 
duction of gravitinos whose decays could aversely affect big bang nucleosynthesis and could 
overpopulate the Universe with dark matter particles [58]. To avoid this overproduction, 
one should require y u < 10 -5 . For the upper limit in this case, we estimate 

y ~ 10" 5 : IV* ~ 51.7, (55) 

as shown by the vertical green line in Fig. [71 On the other hand, as discussed above, the 
smallest value of y consistent with conventional Big Bang nucleosynthesis is y ~ 1CU 16 , in 
which case 

y ~ 10' 16 : IV* ~ 43.4, (56) 

as shown by the vertical purple line in Fig. [71 The above range of couplings includes 
possible gravitational decays of the inflaton |§. We discuss below the compatibility of these 
predictions with the Planck data, shown as the horizontal yellow and blue lines in Fig. | 7 ] 

3.2 Decays via Gravitational-Strength Couplings 

There is another class of no-scale models in which the the inflaton decays via couplings 
that are suppressed by one or more powers of Mp, which we exemplify here by examples in 
which the compactihcation volume modulus T is identified as the inflaton @. For instance, 
in the example discussed in Section 5.2 of [40], there are decays into three-body ttH and 
related final states with rate 

r (T -5 H%t R , i L H%, t R t L H°J = (2 n, + n H - 3) ^,,™ , (57) 

where nt and np are modular weights that are 0(1). Since yt = 0(1) and m ~ 10 -5 Mp, 
this example corresponds to 

3 — body decay : — ~ 5 x 10 -16 . (58) 

m 

§As was discussed in [401 , another possibility in such a matter inflaton scenario would be a superpotential 
coupling of the form £(T — 1/2 ) 2 (j), which would yield decays into T fields with a rate = m|C| 2 /( 367 r). 
The results for different values of y discussed below could also be applied to this case, by simply replacing 
y ->■ Cv^A 

"There are also matter-inflaton models in which decays are suppressed by some power of Mp, but we 
do not discuss them here. 
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In this case we find 


(59) 


3 — body decay : Afi ~ 50.3 , 

as shown by the vertical pale blue line in Fig. |7] 

However, such three-body decays may be dominated by two-body inflaton decays 
into pairs of gauge bosons [261159]. if the gauge kinetic function f a p (where a,/3 are gauge 
indices) has a non-trivial dependence on the volume modulus: f a p = f5 a p with 


d g,T 


(Re/)" 1 



~f~ o, 


(60) 


which is a generic feature of heterotic string effective field theories. In this case, the decays 
into Standard Model gauge bosons V yield 


corresponding to 


r(r ->• vv) 


d 2 T m 3 

y' 1 _ 

32ttMJ, ’ 


Decays into gauge bosons : 



d 2 

u a T 2 

m , 
32t r 


(61) 


(62) 


In a weakly-coupled heterotic string model, one might expect d g T = 0( 1/20), whereas it 
might be 0(1) in a strongly-coupled model, leading to 


Weakly — coupled : — ~ 2 x 10 15 

Strongly — coupled : ^ ~ 10~ 12 . (63) 

These estimates of lead to the following estimates of Ay: 

Weakly — coupled : AT* ~ 50.5 

Strongly — coupled : AT* ~ 51.5, (64) 


as shown by the vertical yellow and magenta lines in Fig. [7] respectively. The compatibility 
of these predictions with the Planck data is also discussed in the next Section. 


4 CMB Bounds on TV* in Representative No-Scale In¬ 
flation Models 

4.1 Matter Inflaton Case 

In the recent no-scale inflation model [29j with an untwisted matter field (j) playing the 
role of the inflaton, the observables (n s ,r) were calculated assuming a Kahler potential 
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of the form (OQ) and choosing a Wess-Zumino superpotential W (0) combining bilinear and 
trilinear terms: 

^ = 2 0 3 03 ’ (65) 

and assuming that the volume modulus T is fixed. It was shown in {29] that this model 
reproduces exactly the predictions of the Starobinsky R + R 2 model if A = /i/3. This model 
can alternatively be written in a more symmetric form: 


K 


— 3 In 



\yi\ 2 + \y2\ 2 ^ ^ 


( 66 ) 


where 


V\ = 


20 


1 + 2T 


2/2 


= Vs 


1-2 T 
1 + 2 T 


in which representation the superpotential (1651) can be written as 


W(y h y 2 ) = l-i 


y\ 


1 + 


2/2 


y\ 


Vs) sVs_ 


(67) 


( 68 ) 


In the Starobinsky limit A = /i/3, and we consider related models with A//i 1/3. 

The calculations of N * in the previous Sections were made assuming exactly Starobinsky- 
like inflation, which (as already mentioned) corresponds in this matter inflation model to 
the limiting case A//i = 1/3. We have studied the modification of the TV* calculation when 
A//r % 1/3, but lying within the range 0.33324 < A//r < 0.33338 displayed in Fig. [8] 

Fig. [H] displays the Planck 2015 constraints on this model in the (n s ,r) plane (upper 
panel) and the (IV*, A//x) plane (lower panel), with the region favoured at the 68% CL 
shaded yellow, and the region allowed at the 95% CL shaded blue. We see in the upper 
panel that for values of A//i ~ 1/3 (black lines) the tensor to scalar ratio is small, and 
in this case, the data yield constraints on n s that are relatively insensitive to r. On the 
other hand, we see that any fixed value of n s corresponds to values of TV* (coloured lines) 
that are strongly correlated with the values of A//i. Thus for a given value of A/p, the 
lower bound on n s provided by Planck can be translated into a lower bound on IV* that is 
sensitive to A//r. For example, for \/y = 1/3, the 68% lower bound on n s corresponds to 
a lower bound of ~ 50 on IV*. 


This feature is reflected in the lower panel of Fig. [HI where we see that IV* is essentially 
unconstrained in this model in the absence of a precise value for A //i. However, if one 
assumes the Starobinsky value A/+ = 1/3, one finds IV* 6 (50,74) at the 68% CL, which 
would disfavour y < 10~ 9 according to Fig. [7] and the 68% CL lower bound on IV* would 
strengthen for A/p > 1/3. 

As seen in Fig. [9] the maximum deviation from the Starobinsky prediction for IV* 
(as shown in Fig. [7]) due to varying A/+ in the range studied (yellow band) is always < 1, 
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Figure 8: The 68% and 95% CL regions (yellow and blue, respectively) in the ( n s ,r ) plane 
(upper panel) and the (IV*, A//i) plane (lower panel) for the no-scale inflationary model 
with a matter inflaton field and the Wess-Zumino superpotential / fi)51) . The black lines in the 
upper panel are contours of X/p, and the coloured lines are contours of N*. The horizontal 
black line in the lower panel is for A//i = 1/3, the value that reproduces the inflationary 
predictions of the Starobinsky model 12$) . 
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and the deviation is significantly smaller for the favoured models with inflaton decay via 
a two-body superpotential coupling y < 10 -5 (corresponding to the green vertical line in 
Fig. [3). The Starobinsky-like analysis gave IV* ~ 51.7 for y = 10~ 5 , as seen in (l55lh and the 
non-Starobinsky deviation of IV* in Fig. [9] is < 0.5 for this value of y , decreasing to much 
smaller values close to the Big Bang nucleosynthesis lower limit y ~ 10 -16 , for which we 
found IV* ~ 43.4 in the Starobinsky limit, as seen in (|56lh 



Figure 9: The possible variation of the value of N* around the prediction for the Starobinsky 
limit X/y = 1/3 in the matter inflation model F6B) with 0.33324 < X/y < 0.33338 (yellow 
band), and for a stabilizing parameter 10~ 2 < A < 1 in the Kahler potential m (blue 
band), as a function of the inflaton decay rate F^. 


It is necessary to address in this model two potential issues: the stabilization of the 
real and imaginary parts of the field y 2 when the inflaton y\ reaches its minimum, and the 
possibility that V,ey 2 0 during inflation. The first of these issues is resolved by adding a 
small supplementary term to the superpotential ()68lk 

A w = bp V i, (69) 

for some constant b. The second issue is addressed by incorporating a quartic term in the 
Kahler potential flfifill : 


K 


— 3 In 


A W + W 



(70) 
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where typical values of A < 1 in natural units, as discussed in [30] • 

Fig. ITOl shows the 68% and 95% CL regions (yellow and blue, respectively) in the (n s , r ) 
plane (upper panel) and in the (AT*, A) plane (lower panel) for this no-scale inflationary 
model with the illustrative choice b = 1CT 6 and the range A < 1. As shown in [30], this 
matter inflaton model reproduces the inflationary predictions of the Starobinsky model 
for A < 1, and we see in the upper panel of Fig. [TO] that the data constraints on n s are 
insensitive to r for A in this range and the relevant values of AT* (coloured lines). In each 
of the segments shown, A varies from 1/100 to 1 as shown for several values of N *. Once 
again, we can use the lower bound on n s to derive a lower bound on Ah for a given value of 
A. Those limits are reproduced in the lower panel of Fig. [TO] where we see that in the limit 
of small A the current data require Ah > 45.5 at the 95% CL and favour Ah £ (51.0, 75.2) 
at the 68% CL. The variation in Ah for values of A £ (10~ 2 ,1) is shown as a thin blue band 
in Fig. [9] as a function of T^/m. We see that this is always smaller than the variation due 
to varying A//i, being <C 1 and hence negligible for our purposes. 

4.2 Volume Modulus Inflaton Cases 

In [38, 39] we calculated the observables ( n s , r ) in various no-scale inflationary models in 
which the inflaton was identified with some combination of the real and imaginary parts of 
the complex volume modulus T. In this case, we choose a superpotential of the form |60j 

W= V3m(j)(T- 1/2). (71) 

Inflation along the direction of the canonically normalized real component p of T yielded 
a Starobinsky-like model, whereas there was a quadratic potential along the imaginary 
direction [MIST]- hi [38] we assumed that higher-order terms in the Kahler potential K 
fixed the angle of the inflationary trajectory in the complex T plane, whereas [39] we made 
a complete two-held analysis of the inhationary observables n s and r. We now confront 
these models with the Planck 2015 data. 

The models are characterized by two parameters, the angle of the starting-point in 
the complex T plane, which we parameterize here as a = arctan(cr/p) 0, and the modulus 
stabilization parameter c [61j: 

K = -3 In (r + T* + c [sin a{T + T* — 1) — cos a(T - T*) 2 ] 2 ) + ^ . (72) 

In this model the matter held </> relaxes dynamically to zero during inflation [38]. As was 
discussed in [38lf39] . the case a = 7t/ 2 corresponds to a quadratic model of chaotic inflation, 

II Note that the angle 6 used in fTSlfoili is equivalent to 7r/2 — a. 
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Figure 10: The 68% and 95% CL regions (yellow and blue, respectively) in the ( n s ,r ) plane 
(upper panel) and the (IV*, A) plane (lower panel) for the no-scale inflatiojiary model ISO |/ 
with a matter inflaton field and the Wess-Zumino superpotential j[68 1 EPj], for b = 10 -6 and 
A < 1. The coloured lines in the upper panel are contours of N* for 0.01 < A < 1. 

and smaller values of a interpolate between this and the Starobinsky limit when a = 0. If 
the quartic stabilization term oc c is large enough, the inflaton trajectory follows a narrow 
valley in held space, much like a bobsleigh run, whereas if c is small the inflaton trajectory 
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is less constrained and two-field effects become important. In this case, we found in [39j 
that the trajectories for a < 7r/2 tend to become more Starobinsky-like for smaller values 
of c. Predictions of this model for various values of c and a, based on a full two-field 
analysis, can be found in [39] . where it can be seen that the results depend on IV*. The 
Planck 2015 68% and 95% CL contours in the (n s ,r) plane can then be converted into 
the corresponding constraints on IV* in these no-scale models as functions of c and a, as 
illustrated in the following figures. 

In Fig. [IT] we display the 68% and 95% CL regions (yellow and blue, respectively) 
in the (n s ,r) plane (upper panel) and in the (IV*, a) plane (lower panel) for the strongly- 
stabilized case c = 100 q We see in the upper panel that in this case the data constraints 
on n s depend in an essential way on the value of r, and we also see that the contours of 
IV* (coloured lines) depend in a non-trivial way on the value of a, as discussed in [39] . The 
interpolation between a = 7t/2, the limit in which the model realizes chaotic quadratic 
inflation, and a = 0, the Starobinsky limit, is nonlinear. 

For these reasons, the conversion of the CMB data into constraints on IV* also depends 
non-trivially on a, as seen in the “whale-like” shape in the lower panel of Fig. [TT] In 
particular, the whale’s “mouth” is the converse of the leftward swerve in the IV* contours 
in the upper panel of Fig. [II] and the “lower jaw” corresponds to the swerve back to larger 
n s at small r. Overall, we see that IV* > 50 is preferred at the 68% CL, whereas IV* > 43 
is allowed at the 95% CL. Translating these limits into constraints on r^/m is also less 
trivial than in the previous matter inflaton models, since this matter inflaton model is not 
Starobinsky-like for large a, as seen in the upper panel of Fig. [TT] As seen in Fig. [12] any 
fixed value of IV* may correspond to a range of values of Tt/iti, shown by the blue band, 
depending on the value of a. The upper side of the band correspond to the limit of chaotic 
inflation with a quadratic potential, and the lower side to the Starobinsky limit. As we see 
in Fig. [7] in the Starobinsky-like limit a —» 0 we find the following constraints on Tt/tti 
and the effective two-body coupling y. 

- > 3 x 10" 20 , y > 10- 9 (68% CL). (73) 

m 

We have studied numerically the possible variations in the dependence of IV* on the decay 
rate Tt in the range 0 < a < 7t/ 2, with the results shown as the blue band in Fig. [12] We 
see there that in the chaotic quadratic inflation case the constraint on Tt/iti is relaxed by 
a factor ~ 10 6 and the effective two-body coupling y is relaxed correspondingly by a factor 
~ 10 3 . 

If y < ICC 5 , corresponding to values of IV* < 52.7 as shown in (|55|) and Fig. [7] for the 
Starobinsky case, then only the range a < 7r/16 is allowed at the 68% CL, rising to a < 7t/4 

**Larger values of c would give very similar results. 
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Figure 11: The 68% and 95% CL regions (yellow and blue, respectively) in the ( n s ,r ) 
plane (upper panel) and the (N*,a) plane (lower panel) for the no-scale inflationary model 
m. assuming strong stabilization with c = 100. The coloured lines in the upper panel are 
contours of IV* for 0 < a < 7t/2 . 


at the 95% CL. Thus models with a starting-point of inflation with small a close to the 
real direction, i.e., close to the Starobinsky model, are favoured in the strongly-stabilized 
case, though not strongly at the 95% CL. In the Starobinsky limit a —> 0, corresponding 
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Figure 12: The variation of the value of 77* as a function of the normalized decay rate 
T r/m for the modulus inflaton model with c = 100 and c = 0 (blue and yellow bands, 
respectively), in the range 0 < a < 7t/2. The stabilized region c = 100 contains fully the 
unstabilized region with c = 0. The upper sides of the bands correspond to the limit of 
chaotic inflation with a quadratic potential, and the lower sides to the Starobinsky limit. 

to the lower jaw of the “whale” in Fig. [lTJ the value y ~ 10~ 5 is at the Planck 2015 68% 
CL limit, whereas IV* > 47 and hence y > 10” 13 is allowed at the 95% CL. 

Fig. m shows analogous results for the unstabilized case c = 0. In this case, as 
seen in the upper panel, the important CMB constraint is that from n s , which is almost 
independent of r for the relevant values of N. * (coloured lines). The lower panel of Fig. fl3l 
shows that again 77* > 50 is favoured at the 68% CL, whereas now IV* > 44 is allowed at the 
95% CL. Both of these limits are rather insensitive to a in the range < 37 t/ 8, whereas only 
very large values of N. * are allowed as a —> tt/2, corresponding to a starting-point along 
the imaginary direction. In this case there is no significant preference for Starobinsky-like 
models with the starting-point of inflation close to the real direction. 

Fig. f T41 shows analogous results for inflation along the imaginary direction, i.e., a = 
7t/ 2, and allowing the stabilization parameter c to vary. As we see in the upper panel, the 
cosmological constraint depends non-trivially on both n s and r, and the values of these 
quantities for fixed N. * (coloured lines) depend non-linearly on c £ [0,1]. In this case, as 
seen in the lower panel of Fig. [T4l the lowest value of IV* is found for c ~ 0.03, with 77* > 48 
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Figure 13: As for Fig. Q71 but assuming no stabilization, i.e., c = 0. 


favoured at the 68% CL and iV* > 43 allowed at the 95% CL. As can be seen in Fig. [I2l 
these constraints correspond to y > 1CF 12 being favoured at the 68% CL, whereas all values 
of y consistent with Big Bang nucleosynthesis are allowed at the 95% CL. This example 
shows that inflation is possible in this model even if the starting-point of inflation is close 
to the imaginary direction, far from the Starobinsky-like limit. 
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Figure 14: The 68% and 95% CL regions (yellow and blue, respectively) in the ( n s ,r ) plane 
(upper panel) and the (IV*, c) plane (lower panel) for the no-scale inflationary model ( f 13) , 
assuming a starting-point along the imaginary field direction. The coloured lines in the 
upper panel are contours of IV* for 0 < c < 1. 


5 Summary and Prospects 


The first purpose of this paper has been to present in Section 2 calculations of N .* in models 
of inflaton decays, with particular attention to predictions in Section 3 from models that 
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are motivated by no-scale models of inflation and have a structure inspired by common 
scenarios for string compactification. We have then analyzed in Section 4 the ranges of IV* 
within no-scale models that yield values of n s and r within the 68 and 95% CL regions 
found by the Planck Collaboration. 

Comparing the results of the previous two Sections, we showed in detail how models 
with smaller inflaton decay rates lead to lower values of T re h, IV* and n s , whereas the data 
prefer larger values of n s and hence IV*. Numerically, in the Starobinsky-like limit models 
with two-body decays and y < 10~ 5 as suggested by the gravitino constraint correspond 
to IV* < 52.7, whereas the data favour TV* > 50 at the 68% CL, with values of IV* > 43 
being allowed at the 95% CL. The present data therefore tend to favour models with 
relatively rapid inflaton decay: T^/m > 1CT 19 corresponding to y > 2 x ICC 9 for two-body 
decays at the 68% CL, to be compared with the upper bound y < 10 -5 from the gravitino 
abundance 0. ft will be interestingto see how, as the experimental constraints tighten, the 
experimental noose on y tightens Ell 

It is interesting to consider specific messages from our analysis for a couple of phe¬ 
nomenological issues, namely sneutrino inflation and supersymmetry breaking. 

The cosmological upper limit on the gravitino abundance imposes an upper limit 
on the Yukawa coupling responsible for sneutrino inflaton decay: y„ < 10 -5 , which is 
itself an important constraint on realizations of sneutrino inflation. If the matter held 
$ in the no-scale Wess-Zumino model (1651) is identified with a sneutrino, it must have a 
trilinear coupling A ~ y/3. Such a trilinear coupling violates A-parity, inducing decay of 
the lightest supersymmetric particle. However, its lifetime is still much longer than the age 
of the Llniverse, and it remains a viable candidate for cold dark matter. As we have shown 
here, if A > /i/3, CMB measurements measurements favour IV* > 50, and hence are on 
the verge of providing a relevant lower bound on the coupling y u responsible for sneutrino 
decay. It will be interesting to see how this squeeze on y u will evolve. 

Concerning supersymmetry breaking, we recall that there is a contribution to gaugino 
masses of the form 


my 2 


^ qG/ 2 fa/3,T 

2 Re f a p 


(G~ l )lG T 


YE|p-3|m3/2 


(74) 


where p is a number of order unity. We saw in Section 3 that strongly-coupled string models 
with d g T = 0(1) would give larger values of the inflaton decay rate T^ ~ d 2 g T m 3 /(32 -k) 
(l62]h T reh , Y* and hence n s than weakly-coupled models with d 9t T = 0(1/20), as seen in 

^However, the lower limit on may be relaxed by a factor ~ 10 6 (and that on y by a factor ~ 10 3 ) in 
non-Starobinsky-like no-scale models. 

Idle mild preference for rapid inflaton decay applies to many models, including the original Starobinsky 
R + R 2 model, no-scale Starobinsky-like models J2S] and Higgs inflation. 
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(EH) . The difference in N* and hence n s is not yet signihcant, but this connection will also 
be interesting to watch in the future. 

These two examples serve as illustrations how inflationary observables may in the 
future join the phenomenological mainstream. These two examples are in the context of 
specific no-scale supergravity models of inflation, but the connection has the potential to 
be more general. 
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